Effective nucleon mass in relativistic mean field theory 
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In the o-uj-p model of the relativistic mean field theory with nonlinear cr-meson self-interaction, 
the effective nucleon mass M* is discussed with relation to the symmetry incompressibility K s of 
nuclear matter, based on the model parameters fitted to nuclear matter properties. It is shown that 
M* is larger than 0.73M, if K s is assumed to be negative and the nuclear matter incompressibility 
Ko is kept less than 300A/eV. Furthermore, the field system is shown to be stable, as the cr-meson 
self-interaction energy is lower bounded in this parameter region. 
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As starting point for the relativistic microscopic description of the nuclear many-body system, within the framework 
of quantum hadrodynamics, the well-studied a-uj-p model with nonlinear cr-meson self-interaction has been proved 
to be able to describe the saturation and other properties of nuclear matter jjj. However, there is some controversy 
CD , about the stability of the field system specified by the fitted parameters, as the fitted coefficient of the fourth power 
term of the cr-meson field in the Lagrangian density is negative, and thus the energy of the field system diverges for 
t-H \ large cr-field [|l| fn addition, the symmetry incompressibility K s of nuclear matter determined by the parameter 
sets existing in the market is positive, which is opposite to that given by the nonrelativistic models of nuclei |}). 
Furthermore, the effective nucleon mass M* is around 0.6M which seems uncomfortably low. Then, the relevant 
question is: are these features manifestation of intrinsic properties of the model itself, or do they depend on the input 
data used for fitting the model parameters? The purpose of this rapid communication is to show that the stable result 
with negative symmetry incompressibility K s and larger effective nucleon mass M* can be obtained, if the model 
■ parameters are fitted to reasonable nuclear matter properties. 

The a-u-p model of the relativistic mean field theory is specified by the following Lagrangian density 0(we use 
natural units with % = c = 1): 
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where — d^io v — d v uo^ , B^" — d^h" — d v h^; -0, <f>, lo and b M are the nucleon, cr, ui and p meson fields with masses 
M, m CT , m w and m p , respectively, while g a , g^ and g p are the respective meson-nucleon coupling constants; b and c 
are the nonlinear term coefficients, and r are isospin matrices. As M, and m p are taken from experiment, the 
model parameters are g a , g u , g p , b, c and m a . 

The nuclear matter equation of state derived from this Lagrangian density can be expressed in terms of the nuclear 
energy density £ as e = £/pn — M, where pn is the baryonic density, and 

£ = £ k + £ a + £u + £ P , (2) 
£ fc = ^ £ iW£M), (3) 
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where 5 = (p p — p n )/pN is the nuclear matter asymmetry; k p and k n are the proton and neutron Fermi momenta 
respectively, 
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and the function F m (x) is defined as (see Rcf. [§ for details): 



F m (x)= dxx 2m ^Jl + x 2 . 



The reduced effective nucleon mass £ and thus the field <p is determined by 

(1 - + bCl{\ - £) 2 + cC 2 a (l - if = 
where the scalar density p s can be expressed as 
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and the function f m {x) is defined as (see Ref. S for details) 

fm(x) = dx 
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The model parameters related to the nuclear equation of state are only C 2 , C 2 , C 2 , b and c. As m a is the inverse 
Compton wavelength of the cr-meson, it is related only to finite size effects, such as surface energy and shell effects of 
nuclei. 

Knowing the equation of state, the following formula for pressure p can be obtained: 



p = -£ + pn^— = - l; M (ps -£a + £u> + £ P - 
opN 3 3 
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The standard state (pn = Pot$ = 0) is defined by the equilibrium condition, p(po,0) = 0. The standard density pg 
can be written in terms of nuclear radius constant rp or nucleon Fermi momentum fcj? as 
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As the finite nuclei are in states near the standard one, the equation of state can be written approximately as 



e(p N , 5) = -ai + -^(K Q + K s 5 2 ) ( 
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which is specified by the standard density po: volume energy ai, symmetry energy J, incompressibility Kq, density 
symmetry L and symmetry incompressibility K s . These 6 quantities can be expressed in terms of 5 model parameters 
C 2 , C 2 , C 2 p , b and c || . Reversebly, the 5 model parameters can be fixed if 5 of those quantities are known. In this 
case, we will choose ro, <zi, Kq, J and K s as input data, where ro is equivalent to po- Specifically, the procedure is as 
follows. 
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At the standard state (po,0), e(po,0) = —a%, and £ p = 0, we have 

£ k + £ a + £ u = po(M-ai). (16) 
In addition, the equilibrium condition, p(po, 0) = 0, can be written as 

^(£ k -M^p s )-£ a + £ u> =0. (17) 
The following formulas can be derived from above two equations: 
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Furthermore, Eqs. (||) and ( |lC)| ) can be combined to solve the parameters & and c as 



& = 1^ ? (20) 
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Finally, the equation for symmetry energy J H is solved to give 
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For given ro, oj, £, C 2 and J, the above formulas can be used to calculate C 2 , C 2 , 6 and c, and thus to calculate Ko, 
K s and L 0. According to the calculated Kg and K s , we can fix £ and C 2 . The following experimentally acceptable 
values |{| are used as input data in our calculation: 

r fa 1.14/m, oi ps l&MeV. (23) 

Fig.l plots the calculated if s versus K for given £. The solid curves from top to bottom correspond to £ = 0.5, 0.6, 
0.7, 0.8 and 0.85, respectively. It can be seen that K s is negative only for £ larger than about 0.7 — 0.8, if Ko < 300 
MeV is assumed. This is shown more clearly in Fig. 2, where the calculated K s versus £ is displayed for given K$. 
The solid curves from top to bottom correspond to Kq — 200, 300, 400 and 500 MeV, respectively. It is shown that 
K s is negative when £ is larger than about 0.73, for K < 300 MeV. 

Fig. 3 is the calculated C 2 versus Ko for given £. The solid curves from top to bottom correspond to £ = 0.5, 0.6, 
0.7, 0.8 and 0.85, respectively. Fig. 4 gives the calculated C 2 versus £. It can be seen from Eq.(|18|) that C 2 depends 
only on £, for given ro and a±. The calculation of C 2 depends, beside £, also on the input value of J and Fig. 5 displays 
the calculated C 2 p versus £ for given ro, ai and J — 30 MeV. On its turn, L depends, beside £ and J, also on the 
input value of K , and Fig. 6 shows the calculated L versus £ for given r , a±, J = 30 MeV and K Q . The solid curves 
from top to bottom correspond to Kq = 200, 300, 400 and 500 MeV, respectively. 

Fig. 7 presents the nonlinear coefficient b versus Ko for given £. On the right hand side of the plot, the first three 
curves correspond to £ = 0.5, 0.6 and 0.7 from top to bottom, respectively. In the middle of the plot, the lower two 
curves correspond to £ = 0.8 and 0.85 from top to bottom, respectively, being the first curve scaled by xl/2 and the 
second one by xl/10. 

Fig. 8 displays the nonlinear coefficient c versus Kq for given £. On the right hand side of the plot, the solid curves 
correspond to £ = 0.5, 0.6, 0.7, 0.8 and 0.85 from bottom to top, respectively. The value of c should be scaled by 
xl/10 for the curve of £ = 0.8, while by xl/50 for the curve of £ = 0.85. It can be seen from Fig. 8 that c is positive 
if £ is larger than about 0.7 - 0.8, for K < 300 MeV. 

In addition to these results, it is worthwhile to see what could be obtained, if the realistic nuclear matter properties, 
extracted from measured data of finite nuclei by nonrelativistic models, are used as input data. In this case, the results 
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given by Myers- Swiatecki phenomenological nucleon-nucleon interaction [|7|, Skyrme interaction M as well as Tondeur 
interaction || are employed. 

The results of the calculation is presented in Table |. The input data set (ro, a\, Kg, J, K s ) is taken from the 
compilation of Ref. 0. MS is for the Myers-Swiatecki interaction, SIII, Ska, SkM, SkM* and RATP for the Skyrme 
interaction, and Tondeur for the Tondeur interaction. It is worthwhile to note that the input value of K s is negative 
for all of these interactions. This is in agreement with the most expectations based on the nonrelativistic model 
p| . As it can be seen from Eq. (|l5|) , physically, K s < means that the compression modulus of nuclear matter 
will decrease if the asymmetry 5 of nuclear matter increases. Experimentally, K s extracted from the isoscalar giant 
monopole resonance energy is between — 566 ± 1350 to 34 ± 159MeV jl^ ]. It is also interesting to note that the effective 
nucleon mass M* /M is around 0.89, which reproduces nicely the value given by nonrelativistic models of nuclei. In 
addition, the nonlinear coefficient c is positive, which means the field system is stable for all of these parameter sets. 

In summary, the effective nucleon mass M* is discussed with relation to the symmetry incompressibility K s of nuclear 
matter, based on the model parameters fitted to nuclear matter properties, in the a-tu-p model of the relativistic mean 
field theory with nonlinear cr-meson self-interaction. It is shown that M* is larger than 0.73M, if K s is assumed to 
be negative and nuclear matter incompressibily Kg is kept less than 300 MeV. Furthermore, it is shown also that the 
field system is stable, as there is a lower limit for the cr-meson self-interaction energy in this parameter region. Our 
conclusion is: a stable result can be obtained by using an appropriate model parameter set altogether with a judicious 
selection of nuclear properties, within the a-u-p model of the relativistic mean field theory with nonlinear cr-meson 
self-interaction, at least for nuclear matter properties. 
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FIG. 1. Calculated K 3 versus Ko for given ro = 1.14 fm, a\ = 16 MeV and £. The solid curves from top to bottom correspond 
to £ = 0.5, 0.6, 0.7, 0.8 and 0.85, respectively. 



FIG. 2. Calculated K 3 versus £ for given ro = 1.14 fm, cii = 16 MeV and Ko. The solid curves from top to bottom correspond 
to K = 200, 300, 400 and 500 MeV, respectively. 



FIG. 3. Calculated C 2 a versus Ko for given ro = 1.14 fm, a± = 16 MeV and £. The solid curves from top to bottom correspond 
to ^ = 0.5, 0.6, 0.7, 0.8 and 0.85, respectively. 



FIG. 4. Calculated C 2 , versus £ for given ro = 1.14 fm and cii = 16 MeV. 



FIG. 5. Calculated C 2 P versus £ for given r = 1.14 fm, en = 16 MeV and J = 30 MeV. 
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FIG. 6. Calculated L versus £ for given r = 1.14 fm, a\ = 16 MeV, J = 30 MeV and K . The solid curves from top to 
bottom correspond to Ko = 200, 300, 400 and 500 MeV, respectively. 



FIG. 7. Nonlinear coefficient b versus Ko for given ro = 1.14 fm, oi = 16 MeV and f . The solid curves correspond to £ = 0.5, 
0.6, 0.7, 0.8 and 0.85, respectively. 



FIG. 8. Nonlinear coefficient c versus K for given r = 1.14 fm, a x = 16 MeV and f. The solid curves correspond to £ = 0.5, 
0.6, 0.7, 0.8 and 0.85, respectively. 



TABLE I. The nuclear matter properties r (fm), oi(MeV), J (MeV), K (MeV), L(MeV), and the parameters C 2 a , Cl, C 2 p , 
b, c of nonlinear a-u)-p model in the relativistic mean field theory. See text for details. 





ro 




Ko 


J 


K s 


L 


M*/M 


cl 


cl 


c P 


b 


c 


MS 


1.140 


16.24 


234.4 


32.65 


-147.1 


85.553 


0.893371 


92.728 


30.908 


27.729 


-0.09203 


1.1137 


SIII 


1.180 


15.86 


355.5 


28.16 


-393.9 


72.866 


0.877399 


77.041 


47.982 


24.665 


-0.15264 


1.0935 


Ska 


1.154 


15.99 


263.1 


32.91 


-78.45 


86.766 


0.885116 


96.522 


38.415 


29.434 


-0.08115 


0.8458 


SkM 


1.142 


15.77 


216.6 


30.75 


-148.8 


79.898 


0.897339 


95.423 


28.962 


25.294 


-0.08287 


1.1499 


SkM* 


1.142 


15.77 


216.6 


30.03 


-155.9 


77.733 


0.897526 


94.984 


28.842 


24.267 


-0.08441 


1.1655 


RATP 


1.143 


16.05 


239.6 


29.26 


-191.3 


75.430 


0.893634 


89.460 


31.269 


23.183 


-0.10318 


1.1808 


Tondeur 


1.145 


15.98 


235.8 


19.89 


-39.78 


47.603 


0.886214 


107.753 


36.378 


9.705 


-0.04911 


0.6885 
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